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Abstract. We consider the Klein-Gordon equation on a star-shaped network composed of n 
half-axes connected at their origins. We add a potential which is constant but different on each 
branch. The corresponding spatial operator is self-adjoint and we state explicit expressions 
for its resolvent and its resolution of the identity in terms of generalized eigenfunctions. This 
leads to a generalized Fourier type inversion formula in terms of an expansion in generalized 
eigenfunctions. Further we prove the surjectivity of the associated transformation, thus showing 
that it is in fact a spectral representation. 

The characteristics of the problem are marked by the non-manifold character of the star- 
shaped domain. Therefore the approach via the Sturm-Liouville theory for systems is not 
well-suited. The considerable effort to construct explicit formulas involving the tunnel effect 
generalized eigenfunctions is justified for example by the perspective to study the influence of 
tunnel effect on the L°°-time decay. 



1. Introduction 

This paper is motivated by the attempt to study the local behavior of waves near a node in 
a network of one-dimensional media having different dispersion properties. This leads to the 
study of a star-shaped network with semi-infinite branches. Results in experimental physics 
[lffi [20] , theoretical physics [15] and functional analysis [H H3] describe new phenomena created 
in this situation by the dynamics of the tunnel effect: the delayed reflection and advanced 
transmission near nodes issuing two branches. It is of major importance for the comprehension 
of the vibrations of networks to understand these phenomena near ramification nodes i.e. nodes 
with at least 3 branches. The associated spectral theory induces a considerable complexity (as 
compared with the case of two branches) which is unraveled in the present paper. 

Perturbation arguments do not seem to be suited for the above mentioned applications, as 
they lead to Neumann series in the solution formulae that hide the interesting effects. 

The dynamical problem can be described as follows: 
Let Ni, . . . , N n be n disjoint copies of (0, +oo) (n > 2). Consider numbers a k ,c k satisfying 
< Cfc, for k = 1, . . . , n and < a\ < a<i < . . . < a n < +oo. Find a vector (u±, . . . , u n ) of 
functions u k '■ [0, +oo) x N k — > C satisfying the Klein-Gordon equations 

[df - c k dl + a k ]u k (t, x) = , k = 1, . . . , n, 

on N\ , . . . , N n coupled at zero by usual Kirchhoff conditions and complemented with initial 
conditions for the functions u k and their derivatives. 
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Reformulating this as an abstract Cauchy problem, one is confronted with the self-adjoint 
operator A = (— • d 2 + ak)k=i ...n m L 2 (N), with a domain that incorporates the Kirchhoff 
transmission conditions at zero. For an exact definition of A, we refer to Section [2j 

Invoking functional calculus for this operator, the solution can be given in terms of 



The refined study of transient phenomena thus requires concrete formulae for the spectral rep- 
resentation of A. The seemingly straightforward idea to view this task as a Sturm-Liouville 
problem for a system (following [27J) is not well-suited, because the resulting expansion for- 
mulae do not take into account the non-manifold character of the star-shaped domain. The 
ansatz used in [27] inhibits the exclusive use of generalized eigenfunctions satisfying the Kirch- 
hoff conditions. This is proved in Theorem 18.11 in the appendix of this paper, which furnishes 
the comparison of the two approaches. 

A first attempt to use well-suited generalized eigenfunctions in the ramified case but without 
tunnel effect [5] leads to a transformation whose inverse formula is different on each branch. 
The desired results for two branches but with tunnel effect are implicitly included in [27] . For n 
branches but with the same and a& on all branches a variant of the above problem has been 
treated in [7J using Laplace transform in t. 

In the present paper we start by following the lines of [S]. In Section El we define n families 
of generalized eigenfunctions of A, i.e. formal solutions Fj? for As [ai, +oo) of the equation 



satisfying the Kirchhoff conditions in zero, such that e ±ty ^ t F^(x) represent incoming or outgoing 
plane waves on all branches except Nk for A £ [a n , +oo). For A 6 [a p , a p +i), 1 < p < n we have 
no propagation but exponential decay in n — p branches: this expresses what we call the multiple 
tunnel effect, which is new with respect to [5]. Using variation of constants, we derive a formula 
for the kernel of the resolvent of A in terms of the F£. 

Following the classical procedure, in Section H] we derive a limiting absorption principle for 
A, and then we insert A in Stone's formula to obtain a representation of the resolution of the 
identity of A in terms of the generalized eigenfunctions. 

The aim of the paper, attained in Section is the analysis of the Fourier type transformation 



in view of constructing its inverse. We show that it diagonalizes A and determine a metric 
setting in which it is an isometry. This permits to express regularity and compatibility of / in 
terms of decay of Vf. 

A major task in this context is to overcome the cyclic structure that the cyclic nature of 
the n-star induces in the underlying resolvent formula derived in Section [3j To this end, we 
use a symmetrization procedure that is carried out in Section [5j It combines the expression 
for the resolution of the identity E{a, b) found in Section 0] with an ansatz for an expansion in 
generalized eigenfunctions: 



This creates a (3n 2 + 1) x n 2 linear system for the qi m , whose solution leads to the result in 
Theorem 15.31 and to the Plancherel type formula. 




uq and e 
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A direct approach to the same symmetrization problem, carried out in Section [6l yields a 
closed formula for the matrix q based on n x n matrices. This approach is to a great extent 
independent of the special setting and is thus supposed to be generalizable. 

In Section [7] the desired inversion formula as well as the Plancherel type theorem are stated. 
Finally the domains of the powers of A are characterized using the decay properties of Vf. We 
show that V is an ordered spectral representation (see Definition XII. 3. 15, p. 1216 of |17|). The 
spectrum has n layers and it is p-fold on the frequency band [a p ,a p +i). This reflects a kind 
of continuous Zeemann effect caused by the constant, semi infinite potentials on the branches 
Nj given by the terms cijUj. On this frequency band the generalized eigenfunctions have an 
exponential decay on n — p branches, expressing the multiple tunnel effect. 

Our results are designed to serve as tools in applications concerning the dynamics of the tunnel 
effect at ramification nodes. In particular, we think of retarded reflection (following [8j 120] ). 
advanced transmission at barriers (following [18l [151 E] ) , £°°-time decay (following [3]), the 
study of more general networks of wave guides (for example microwave networks [25] ) , causality 
and global existence for nonlinear hyperbolic equations (following [9j) and the generalization to 
coupled transmission conditions (following [12]). 

Finally, let us comment on some related results. The existing general literature on expansions 
in generalized eigenfunctions 1241 127] for example) does not seem to be helpful for our kind 
of problem: their constructions start from an abstractly given spectral representation. But in 
concrete cases you do not have an explicit formula for it at the beginning. 

In [1U] the relation of the eigenvalues of the Laplacian in an L°°-setting on infinite, locally finite 
networks to the adjacency operator of the network is studied. The question of the completeness 
of the corresponding eigenfunctions, viewed as generalized eigenfunctions in an L 2 -setting, seems 
to be open. 

In [22], the authors consider general networks with semi-infinite ends. They give a construction 
to compute generalized eigenfunctions from the coefficients of the transmission conditions and 
the asymptotic behaviour of eigenvalues is studied. The generality of the aproach, however, does 
not allow for explicit inversion formulas for a given family of generalized eigenfunctions. 

Spectral theory for the Laplacian on finite networks has been studied since the 1980ies for 
example by J. P. Roth, J. v. Below, S. Nicaise, F. Ali Mehmeti. A list of references can be found 
in p. 

In [23] the transport operator is considered on finite networks. The connection between the 
spectrum of the adjacency matrix of the network and the (discrete) spectrum of the transport 
operator is established. A generalization to infinite networks is contained in [16J. 

For surveys on results on networks and multistructures, cf. [H 119]. 

Many results have been obtained in spectral theory for elliptic operators on various types of 
unbounded structures for example [211 H31 13 [3] , cf. especially the references mentioned in [3]. 

2. Data and functional analytic framework 
Let us introduce some notation which will be used throughout the rest of the paper. 

Domain and functions. Let N\, . . . , N n be n disjoint sets identified with (0, +oo) (n € N, 
n > 2) and put N := U&=i Nk, identifying the endpoints 0, see [T] for a detailed definition. 
Furthermore, we write [a, &]jv fc for the interval [a, b] in the branch For the notation of 

functions two viewpoints are used: 

• functions / on the object N and fk is the restriction of / to Nk- 

• n-tuples of functions on the branches N}~; then sometimes we write / = (fi, . . . , f n ). 
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Transmission conditions. 

« 

(T ): (u k ) k =i,...,n G C°(N k ) satisfies it»(0) = u fc (0) V i, & G {1, . . . , n}. 
fc=i 

This condition in particular implies that {u k ) k= % ... n niay be viewed as a well-defined function 
on iV. 

n n 

(Ti): ('Ufc)fe=i J ...,„ € Y\ ^(Nk) satisfies ^c fc • d x u k (0 + ) = 0. 
fc=i fc=i 

Definition of the operator. Define the real Hilbert space H = L 2 (N k ) with scalar 

product 

(«,«)# = ^(^,«fc)L2(AT fc ) 

k=l 

and the operator A : D(A) — > H by 

n 

D(A) = |(itjfc)fe=i,... )n G JJiI 2 (iV fc ) : («fe)jfe=i,...,n satisfies (T ) and (Ti)|, 

fc=i 

^((■Wfc)fc=l,...,n) = {A k Uk)k=l,...,n = (-Cfc • dlu k + a k Uk)k=l,...,n- 

Note that, if c k = 1 and = for every k G {1, . . . ,n}, A is the Laplacian in the sense of the 
existing literature, cf. pTQ], [22] . 

Proposition 2.1. The operator A : D{A) — > H defined above is self-adjoint and satisfies o~(A) C 
[at, +oo). 

Proof. Consider the Hilbert space 

V = < {u k ) k =x,...,n G \\ ^(Nk) : (uk)k=i,...,n satisfies (T ) \ 
^ k=i ' 

with the canonical scalar product (■, -)y. Then the bilinear form associated with A + (e — a\)I 
is a £ : V x V — > C with 

n 

a e (u,v) = y~] [c k (d x u k , d x v k ) L 2( Nk ) + (a k + e- ai){u k ,v k ) L 2^ Nk )]. 

k=l 

Then clearly there is a C > with a e (u,u) > C(u, u)y for all u G and all e > 0. By partial 
integration one shows that the Friedrichs extension of (o e , V, .ff) is (^4 + (e — ai)/, -D(A)). Thus 
the operator A + (e — cti)/ is self-adjoint and positive. Hence a(A + (e — ai)I) C [0, +oo) for all 
e > 0, what implies the assertion on the spectrum. □ 



3. Expansion in generalized eigenfunctions 

The aim of this section is to find an explicit expression for the kernel of the resolvent of the 
operator A on the star-shaped network defined in the previous section. 

Definition 3.1. An element / G Y\ k= \ C°° (N k ) is called generalized eigenfunction of A, if it 
satisfies (To), (Ti) and formally the differential equation Af = Xf for some A G C. 
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Lemma 3.2 (Green's formula on the star-shaped network). Denote by V\ t ... i n the subset of the 
network N defined by 

n 

V h ,...,i n = {0}u(j(0,l k ) Nk . 

k=l 

Then u, v G D(A) implies 



n n 

u"{x)v{x)dx= \ u(x)v"(x) dx — ^~]u(lk)v'(lk) + u'(lk)v[lk)- 

'„ JVu 1~ 1 1 1 1 



V h,-,ln J Vl 1 ,...,l n k=l k=l 

Proof. Two successive integrations by parts are used and since both u and v belong to D(A), 
they both satisfy the transmission conditions (To) and (T\). So 

n n 

X>*(oK(o) = «i(o)$>Uo) = o. 

k=l k=l 

Idem for ELl U U°)«*( )- D 

This Green formula yields now as usual an expression for the resolvent of A in terms of the 
generalized eigenfunctions. 

Proposition 3.3. Let A G p(A) be fixed and let e\, e\ be generalized eigenfunctions of A, such 
that the Wronskian tt^ 2 (x) satisfies for every x in N 

w$ >2 (x) = det W(£{x),4{x)) = e\{x) ■ (e^)'(x) - (e^)'(x) • e* (x) ? 0. 

If for some k G {1, . . . , n} we have e^jv^ S H 2 (Nk) and e\\N m G H 2 {N m ) for all m ^ k, then 
for every f G H and x G Nk 



[R(X,A)f](x) 



e\(x)el(x')f{x') dx' + / el{x)e\ (x')f(x') dx' 

(x,+oo) Nk J N\{x,+oo) Nk 



Cfe^l j2 (^) 

Note that by the integral over N, we mean the sum of the integrals over Nk, k = 1, . . . , n. 

Proof. Let A G p(A). We shall show that the integral operator defined by the right-hand side of 
(PQ) is a left inverse of XI — A. Let u G D(A) and x G Nk- Then 



/ A := / e^(x)e^(x , )(A/ - ^)u(x') dx' + / e$(x)e\{x')(\I - A)u{x') dx' 

J(x,+oo) Nk JN\(x,+oo) Nk 

= e\(x) lim f k el{x'){XI - A)u{x) dx + e$(x) lim / e^(x)(AJ - ^)n(x') fix, 

Jx l m ^oo,m^k J Vh lk _ liXtlk+1 ,...j n 

due to the dominated convergence Theorem, the integrands being in L X (R) by the hypotheses. 
We have u G D{A) C \[ n j=l H 2 (Nj) and 

^2 1 iVfc G ^ 2 (iVfc), e^|Ar m G # 2 (iV m ), m ^ & 

by hypothesis and thus 



dstiliVfcOc) • e$\ Nh (x) ^-^-^0, u\ Nk (x) ■ d x e 2 \ Nk (x) ^^^0, 

d x u\ Nm {x) ■ e\\ Nm (x) —> 0, u\ Nm (x) ■ d x e\\ Nm (x) — > 0,m^k, 

x— >+oo x— >+oo 



all products being in some H (Nj). Recall that 

' /"« = / 6 /</' - /(%'(&) + /'(&)</(&) + /(*)</ to - /»</(<0 
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for f,g G H 2 ((a, b)). Now Lemma [3721 and (XI - A)e^ = for r = 1, 2 imply 
I x = e\(x) lim [ I k [XI - A)el(x')u(x') dx' 

/fc^OO LJ X 

+ c fc (-u(l k )d x e%(l k ) + d x u(lk)e%(l k ) + u(x)d x .e2 (x) - d x u(x)e2 (x) 



(x)[ lim / (XI - A)e\(x')u(x') dx 1 

u v l 1 ,-..,l k _ 1 ,x,l k+1 ,...,l„ 



+ y~]cj lim ^-u^d^^) + 0a.it(Z 3 -)e^(Zj)j + c fc ^-u(x)5 !r e^(x) + d x u(x)e\(x) 

= cJe^x^e^x) - d x e\(x)e%(x)\u(x) 
= c k w^ 2 (x)u(x). 

Now the invertibility of XI — A implies the result. □ 
Definition 3.4 (Generalized eigenfunctions of ^4). For k G {1, . . . ,n} and A G C let 



t,k(X) ■= \ and s k := r — . 



Here, and in all what follows, the complex square root is chosen in such a way that V r ■ e 1 ^ = 
w i t h r > o and cf> G [— 7r,7r). 
For A G C and j, k G {1, . . . , n}, F^ j : N -> C is defined for x G 1% by J '(x) := i^jjf (a?) 
with 

J F \,j( x ) = cos(^(A)x) ± fs j (A)sin(^-(A)x), 

I F \,k( x ) = exp(±i^(A)x), for ^ j. 

Remark 3.5. • -Fjf' 5 ' satisfies the transmission conditions (To) and (Ti) and formally it 

holds AF X = XF^. Thus it is a generalized eigenfunction of A, but clearly F x J does 
not belong to H, so it is not a classical eigenfunction. 
• For Im(A) 7^ 0, the function F^£, where the +-sign (respectively — sign) is chosen if 
Im(A) > (respectively Im(A) < 0), belongs to H 2 (N k ) for k 7^ j. This feature is used 
in the formula for the resolvent of A. 

Definition 3.6 (Kernel of the resolvent). Let w : C — > C be defined by w(X) := ±i-^™ =1 Cj^j(X). 
For any A G C such that w(X) 7^ 0, j G {1, . . . , n}, and for every x G Nj, we define 



K(x, x , A) 



—F±f(x)F±'P +1 (x'), for x' G iV,-, x' > x, 



u,(A) 



1 

-F X j + (x)F x ,3 (x f )i for x G iVfc, 7^ j or x' E iVj, a/ < x. 



I ™(A)' 

In the whole formula + (respectively — ) is chosen, if Im(A) > (respectively Im(A) < 0). Finally, 
the index j is to be understood modulo n, that is to say, if j = n, then j + 1 = 1. 

Figure Q] shows the domain of the kernel K(-,-,X) in the case n = 3 with its three main 
diagonals, where the kernel is not smooth. We will show in Theorem 13.81 that K is indeed 
the kernel of the resolvent of A. In order to do so, we collect some useful observations in the 
following lemma. 
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Figure 1. NxN in the case n = 3 

Note that in particular, if Cj = c and a,j = for all j G {1, . . . , n}, then w(X) = ±otc\/A for 
all j € {1, . . . ,n}, which only vanishes for A = 0. On the other hand, if there exist i and j in 
{1, . . . , n}, such that aj ^ a,j, then it is clear that w(X) never vanishes on R, but we need to 
know, if it vanishes on C. 

Lemma 3.7. i) For a\ < X and e > holds \w(X — ie)\ 2 > Y^Jj=i c j\^ ~ a j\- 

ii) For X G p{A) such that Re(A) > a%, the Wronskian w only vanishes at A = a, if = a 
for all k G {1, . . . , n}. 

Proof. We first prove E]). Note that for z\, Z2, ■ ■ ■ , z n G C holds 

j=l j=l k,l=l 



With Zk := Ck£,k{X — ie) and the abbreviation rjk := \JX — ie — ak it follows 



\w{X - ie)\ 2 = ^cjlr/jl 2 + 2 CkCiRe(n k r]i)- 

7=1 fe,!=l 

Thus it suffices to show Re(%r/[) > for /c,Z = l,...,n with k ^ I. With our convention 

this means 

axg(rj k fji) G 



, - arg(z) 

y\z\e l 2 ; arg(z) G [— vr,7r), this means 



7T 7T~ 

2' 2. 

Without loss of generality, let k < I. Then we have < a/ and there are three possible positions 
of A: 
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ai < A: 

Then for r = k and for r = I we have arg(A — ie — a r ) G [— ^, 0] and therefore 



arg(?7 r ) = ^ arg(A - ie - a r ) G 
Using A — afc > A — ai, this implies 

&Tg(j] k r}i) = avg(r] k ) - arg(rft) G 



7T 



.0 



7T 



• cik < A < a^: 
Then arg(A — ie 



afc) G [— f,0] and therefore 
arg(%) = - arg(A - ie - a fc ) G 



Furthermore, arg(A — ie — aj) G [— tt, — § ], and therefore 

1 



arg(r?z) = ^ arg(A - ie - a{) G 
Putting everything together we get 

argimm) = arg(%) - arg(^) G 



7T 7T 

2'~4 



0. 



7T 



A < a k : 

In this case we get, again for r = k and r = I, arg(A — ie — a r ) G [— tt, — ^] and therefore 



1 



arg(r? r ) = - arg(A - ie - a r ) G 



7T 7T 

2'~4 



7T 7T 

T 4 



This yields 

arg(%r?[) = arg(r/ fc ) - arg(^) G 

Thus, in all three cases we have arg^?^) G [— f , f] and, hence, Re^r^) > for all k, I = 
1, . . . , n with k ^ I. 

In order to prove |n]) , we note that the choice of the branch cut of the complex square root 
has been made in such a way that yfpZ = ^fpZ for all p G C. Thus w(~p) = w(p) for all p G C. 
This implies that the first part of the lemma can be generalized to: 

n 

\w(p)\ 2 > ^2cj\Re(n) - aj\ 
i=i 

for every p such that Re(/i) > a\. Then w(p) = for such a /i implies |Re(/i) — aj\ = for 
all j G {1, . . . , n}. Thus all a&, A; = 1, . . . , n have to be equal to some number a > 0. Then 
|u;(^t)| 2 = (Ylk=i y/ck) If 1 ~ a li which only vanishes for // = a. □ 
Theorem 3.8. Let f £ H. Then, for x G N and A G p{A) such that Re(A) > a\ 

[R(X,A)f](x)= [ K(x,x',X)f(x')dx', 

Proof. In (pQ), the generalized eigenfunction e\ can be chosen to be F\'^ ■ Then e\ can be F^' 1 
with any I ^ j so we have chosen j + 1 to fix the formula. The choice has been done so that 
the integrands lie in L 1 (0, +oo) (cf. the last item in Remark 13, 5p . Simple calculations yield the 
expression for the Wronskian w(X). 

For A G p(A) such that Re(A) > ai, the Wronskian only vanishes at A = a if a k = a for all 
k G {1, . . . , n} due to Lemma [3771 But in this case, the Wronskian is w(X) = (Ylk=i yf^k) VX — a 
and u>(A) -1 has an L 1 -singularity at A = a. □ 
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4. Application of Stone's formula and limiting absorption principle 
Let us first recall Stone's formula (see Theorem XII. 2. 11 in |17j). 

Theorem 4.1. Let E be the resolution of the identity of a linear unbounded self-adjoint operator 
T : D(T) -> H in a Hilbert space H (i.e. E(a,b) = l( 0)b) (T) for (a,b) E M. 2 , a < b). Then, in 
the strong operator topology 

h(T)E(a, b) = lim lim — / h(X) [R(X - ei, T) - R(X + ei, T)] dX 
s^o+ e^o+ 2m J a+5 

for all (a, b) E R 2 with a < b and for any continuous scalar function h defined on the real line. 

To apply this formula we need to study the behaviour of the resolvent R(X, A) for A approach- 
ing the spectrum of A. Lemma 13.71 will be useful as well as the following results. First we state 
an estimate that follows directly from the definition of Sj. 

Lemma 4.2. Let 5 > be fixed. For all a± < X, < e < 5 and j = 1, . . . , n holds 

1 - 

\ Sj (X-ie)\<M(X,5):= max { _ . ) V((A - a k f + 5 2 f'\ 

Note that M(-, 5) E L\ oc {[ai, +oo)). Furthermore, if a\ = . . . = a n , then 

1 n _ 

V c i fc=i 

for all fj, E C, which means that Sj is constant and we may take 

n 



M(X,S):= max {— ^^1. 

j=l,...,n I ^fcj j-i J 



Theorem 4.3 (Limiting absorption principle for A). Let 5 > be fixed and let M(X,5) be 
defined as in Lemma \4.3\ Then for all a± < A, < e < <5 and (x,x') E N 2 we have 

i) liniQ.-s.o K(x, x' , X — ia) = K(x, x' , X), 

ii) \K{x,x',X-ie)\ < N{X,5)e^ x+x '\ where N{X,5) := _ n 1+A f| A ' ,5) l , 1/2 and 

c j\*~ a j\) 

_ i x 
7 := max J= i r .. ! „,{c J 2 }max{((a n - ai) 2 + 5 2 )i, 1,5}. 

Proof. i) The complex square root is, by definition, continuous on {z E C : Im(z) < 0} (cf. 

Definition I3.4p . hence the continuity of K(x, x' , •) from below on the real axis. Note that 
x and x' are fixed parameters in this context, 
ii) For fi = X — ie and x E Nj we have in concrete terms 

{[cos(£,j([i)x) — isj(^) sin(£j(/i)x)] exp(— i^j(ii)x'), x E Nj, x > x, 
exp(— i£j(n)x) [cos(^j(/x)j; / ) — isj(^) sin(^(^)x')] , x E Nj, x < x, 
exp(— i^j(n)x) exp(-i£ fc (/i)a/), x E N k , j. 

Now, let us first look at the case A > a n . Then 

\exp(-i^j(fj,)x)\ < exp(|Im(^(/i)x)|) = exp(cj 1/2 |lm( y/X - ie - aj)\). 

Using the fact, that for z E C with | arg(z)| < tt/2 we have lln^-y/i)! < max{l, |Im(z)|}, 
we obtain 

| exp(— i^j(fj,)x)\ < exp(cj 1//2 max{l, 5}x) . 
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In the case a\ < A < a n we find 



\W\ = J^=^ = cT 1/2 ((A - a,) 2 + e 2 ) 1/4 < cT^ ((an _ aj f + ^ 



Using these estimates and |e z | = e Re ^-* < e' 2 ', we find for all A > a\ 

(1 + \8j{ji)\) exp^O^expfl^)^'), x' G iVy, 
exp(|£j(//)|a;)exp(|f fc (//)|a;'), a;' G N k ,k ^ j 




\K(x,x',fj)\ < 

|'uy W | 

- ^^T( 1 + l s j(/ i )l) ex P(7(^ + ^))- 

Hm)I 

The conclusion now follows using Lemma 13.71 and Lemma 14.21 □ 

Note that these estimates in particular imply that N(-,5) G Lj^Qai, +oo)). In fact, if a\ = 
. . . = a n , then M(X,5) can be chosen to be constant, see Lemma 14.21 and the denominator of 
N causes only an Lj oc -type singularity. On the other hand, if there are two different a,j, the 
denominator of N is never zero and M(-, 5) G Lf ([ai, +oo)) again by Lemma 14.21 



Lemma 4.4. For (x,x') G N 2 and X G C, it holds K(x,x',X) = K(x,x',X). 

Proof. The choice of the branch cut of the complex square root has been made such that Va = 
VX for all A G C. This implies e^ x = e iy/lx = e -^ x for all A G C and x G M. Thus it holds 

F+' j (x) = F J ' j (x) and F~' j (x) = F+' j (x) 



for all A G C, x G N and j G {1, . . . , n}. In the same way we have w(X) = —w(X). Observe, 
that switching from A to A the sign of the imaginary part is changing, so in the definition of 
K(x, x' , A) we have to take the other sign, whenever there is a ±-sign in the formula. This gives 
the assertion. □ 

Now, we can deduce a first formula for the resolution of the identity of A. 

Proposition 4.5. Take f G H = Tlj=i L 2 (Nj), vanishing almost everywhere outside a compact 
set B C N and let — oo < a < b < +oo. Then for any continuous scalar function h defined on 
the real line and for all x G N 

(h(A)E(a,b)f)(x) 



L 



h(X) £ lj^(x){ / f(x') l ( ,, +oo) OO • hn(-—F-'\x)F^ +1 (x') 

(a,b)n oi,+oo) ~i J l J 1 3 \W{A) 

+ IaAOe.+oo)^ ix') ■ Im^-^yF A _ ' J+1 (x)F A _J (x / )) dx'} dX, 

where E is the resolution of the identity of A (cf. Theorem \4-l\ ) and the index j is to be 
understood modulo n, that is to say, if j = n, then j + 1 = 1. 

Proof. The proof is analogous to that of Lemma 3.13 in [2]. Let g G H be vanishing outside B. 
Then 

(h(A)E(a, b)f, g) H = ( lim lim -L f & h (X) [R(X - ei, A) - R(X + ei, A)] dX /, g] (2) 
V<5^o+ e^o+ 2m J a+S J h 

= lim lim — ( f h(X)\R(X-£i,A)-R(X + si,A)] dXf,g] (3) 
5^0+ e^o+ 2tti \J a+s J H 



THE MULTIPLE TUNNEL EFFECT ON A STAR-SHAPED NETWORK 11 

lim lim — / h(X) ( \R(X - ei, A) - R(X + ei, A)] /, g) „ dX (4) 
e^o+ 2m J a+S 

hWyJ f(x')[K(-,x',X-ie) - K(-,x',X + ie)] dx',g(-)^J dX 



lim 

<5,e->0+ 2m 



(a+5,fe-(5)n[ai,+oo) 



(5) 

1 



: lim — / h(X)( / f(x') [K(-,x, X — ie) — K(-, x', X - iej] dx', g(-) dX 

8,e^0+ 2m J ( a+ S,b-5)n[ai,+oo) \J N /H 

(6) 

: lim ^- [ h(X)( I f(x')2ilm(K(;x',X-ie)) dx',g(-) ] dX (7) 

S,e->0+ 2m J ( a +8,b-5)n[ai,+oo) \J N Jh 

lim - [ h(X)( I f(x')[ lim Im(K (■ , x' , X - ie))} dx',g(-) j dX (8) 

\J N e-+0+ J H 



8^0+ vr 
/ 1 



(a+<5,&-5)n[ai,+oo) 



7T 



(a,6)n[ai,+oo) JiV 



h(X) f(x')lm(K(-,x' ,X-iQ)) dx' dX,g(-)\ (9) 



/ 1 / h(X){[ /(^[^^^^(l^^^XxOF-^x)^-^ 1 ^) (10) 

7iV 7T 7( a ,6)n[ai,+oo) W7V L ™(A) f-f J V )N 3 



w(X) 4-? 

.7 = 1 



+ l7V\( a; ,+oo^.(^)^ J+1 W^' J ( a;, ))] ^A<?(x) dX. 

Here, the justifications for the equalities are the following: 
([2|): Stone's formula (Theorem I4.ip. 

([3|): After applying the operator valued integral to /, the two limits are in H. So they 
commute with the scalar product in H. 

{'fid)H ls a continuous linear form on C(H), and can therefore be commuted with 
the vector-valued integration. Note that A i— > R(X, A) is continuous on the half-plane 
{A G C : Re(A) < ai}, since the resolvent is holomorphic outside the spectrum, cf. 
Proposition 12.11 

©: Theorem 1^81 

©: Lemma Ol 

(0): z-z = 2i- Imz Mz G C. 

(J8j) : Dominated convergence. Since supp /, supp g and [a, b] are compact, we use the limiting 

absorption principle (Theorem I4.3j) . 
([9]): Fubini's Theorem. 
(UJJ]): Definition ESI 

The assertion follows, since g was arbitrary with compact support. □ 

The unpleasant point about the formula in the above proposition is the apparent cut along 
the diagonal {x = x'} expressed by the characteristic functions in the variable x' . In fact, there 
is no discontinuity and the two integrals recombine with respect to x' . This is a consequence of 
the next lemma that gives an explicit representation of the integrand above. 

In the following, we use the convention 

a n +i ■= +oo, 

in order to unify notation and we set 
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Lemma 4.6. Let j,k G {1, . . . , n} and let A be fixed in (a p , a p+ \), with p G {1, . . . , n}. Then 
Im[^(F A ( X )(F\ J ')fc( a; 0] * s ff^en by the following expressions, respectively: 

• Ifk>j>porj>k>p (Case (a)) 

• Ifj<k<pork<j<p (Case (b), j / fc) 

Im (u;) cos ^' x ) cos (£kx') - Im (^) sin(^x) sin(£ fc x') 

— Re^— ^ cos(£j#) sin(£fcx') — Re^— ^ sin(£jx) cos(£fcx'), 

• If 3 = k<p (Case (b), j = fc) 

lm iw) cos (^' x ) cos (6»') - Im (~) sin (tj x ) sin fe^0 

— Ref— cos(^x) sin(^fcX / ) — Ref — ") sin(£jx) cos(^x'), 



• // j <p<k (Case (c)) 

Im^-je"^ x ' cos(£,-x) + Im^je~^ x ' sin(^-x), 

• Ifk<p<j (Case (d)) 

Imf— )e~^ x cos(£ k x') - lm( — ) e~^ x s in(£ fe x'). 
ViD/ \iwJ 

Proof. Since A belongs to (a p ,a p+ i), £j(A) is a real number, if and only if j < p. Otherwise, £j 

is a purely imaginary number and we have £j = 

Now, the proof is pure calculation, using the following expressions for the generalized eigen- 

functions in the case j / k 

Re(F-'') fe (x) = hff^ if ^ G M ' Le - k ~ P > 



Im(F-^) fe (x) 



e *fc x , if £fc G z'R, i.e. k > p, 

-sin(£ fc x), if £ fc G R, i-e. & < p, 
0, if ^ G zR, i.e. A; > p, 



and for j = k 

R fp-^ C ^ - / cos fe x ) + Im ( s fc) sin (^ x )' if^GR, i.e. fc<p, 

ne ^A J fc W | Re (i (1 + Sfc) ) e -€^ + Re _ flfc ^ ^ fc G iR) j e fc > p) 

, x f -Re(s fc )sin(^x), if £ fc G R, i.e. k < p, 

1 A \ Im (i(l + Sfc )) e"^ + Im (1(1 - s fc )) e** x , if £ fc G iR, i.e. k > p, 

respectively. □ 

Theorem 4.7. Take f £ H = YYj=i ^ 2 (-^j)> vanishing almost everywhere outside a compact set 
B C N and let — oo < a < b < +oo. Then for any continuous scalar function h defined on the 
real line and for all x G N 

(h(A)E(a,b)f)(x)= [ h(X)J2^){f f(x')Im\^F^ + \x)F^(x>)}dx'} d\, 

J[a,b]n[ai,+oo) ~[ 1 l JN lw{\) J J 

(11) 

where again the index j is to be understood modulo n, i.e, if j = n, then j + 1 = 1. 
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Proof. All the work has already been done. It remains to inspect the formulae for j = k in the 
cases (a) and (b) of Lemma 14.61 to observe that the expressions are symmetric in x and x'. So 
for x,x' G Nj with x < x' we find F^' J 1 (x)F7^ +l [xf) = F x '\x')F x which implies the 

assertion. □ 

5. Symmetrization 

As was already explained in the introduction, the aim of this section will be to find complex 
numbers qi m , l,m € {1, . . . , n}, such that the resolution of identity of A can be written as 

(E(a,b)f) (x) = f QimW^x) / F^(x f )f(x f ) dx> dX, (12) 
Ja l,m=l JN 
in order to eliminate the cyclic structure of the formula in Proposition 14.71 

In this section we shall often suppress the dependence on A of several quantities for the ease 
of notation, so sj = Sj(X), q l>m = qi, m (X), fj = £j(A), = ^(A) and w = w(X). 

Lemma 5.1. Given x G Nj, equation (|12 j) is satisfied for all ai < a < b < +oo and all 
f € L 2 (N) with compact support, if and only if for all j = 1, . . . , n 



Im 



-F 



-j+i, 



x)F- 



w 



q lm (X)F-' l (x)F-' m (x') 



(13) 



Z, 771=1 



for almost all x' £ N and X > ai . Here again the index j has to be understood modulo n, that 
is to say, if j = n, then j + 1 = 1. 



Proof. If functions qi t m, I, m = 1, . . . , n, satisfy (fT3|) and if a\ < a < b < +oo, we get by (fTTj) 

\ n „ n 

(E(a,b)f)(x) = j ^l^aOj j fix') ^ q lm {X)F^ l {x)F^ m {x') dx'} dX 

a i= l N l > m=1 

= T E <lim{X)F-\x) [ F~^(x')f(x')dx'dX, 
Ja l,m=l JN 

which is (fT2j) . 

For the converse implication, let (|12p be satisfied for some x G iV,- and all a\ < a < b < +oo, 
as well as all / G L 2 (N) with compact support. This means by (jlip 




a JiV 



/(x )Im 



1 



F7^\x)F^{x') 



Vw{X) 



dx'dX 



b n 



l,m=l 



Firstly, we want to see that the integrands of the A-integrals on both sides in fact have to be 
equal almost everywhere. In order to do so, we observe that they both are in L 1 ((ai, +oo)) and 
use the following general observation: If I is an interval and g G L l (I) satisfies fjg = for 
all intervals J C J, then g = almost everywhere in I. Indeed, in this case we have for any 
Lebesgue point xq G I of g and every e > (cf. [261 Theorem 8.8]) 



\9(x )\ 



^ rx +£ 



< 



2t 
1 

2s 



g{xo) - g{x)) dx + 



X -£ 



2 £ : 



g(x) dx 



XQ+E 



\g(x) - g(x )\ dx — ^0 (e->-0), 



X -£ 
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which implies g(xo) = for almost all xo G /. 
This implies 



N 



f(x )Im 



1 



w(X) 



F^+\x)F^(x') 



dx 



N 



£ qim WF-\x)F^ m {x')f{x')dx' 



Lm=l 



for almost all A > ai and all / G L 2 (N) with compact support. By the fundamental theorem of 
variational calculus this implies the assertion. □ 

In a next step, we explicitely write down equation (fT3|) as a linear system for the values qi m . 

Lemma 5.2. The equation 



Im 



^(F-^).(x)(F^) k (x') 



£ <UA)(F-'') (*)(F->™) fc (x') 



i,m=l 



ZioWs /or any G iVj x AT fc and A G (a p ,a p _|_i), i/ and only if 

• Case (a): if k > j > p or if j > k > p 

Qjk = 

Ez#<fifc = 

Em^A; Qjm = 
J2l^j,m^k llm = Im [-] , 

Case (b), j ^ k: ifj<k<porifk<j<p 



Y.l^j,m^k llm +Sk HljLj Ilk + Sj Ylm^k Qjr 

J2lytj,m^k llm +Sk E/^j Ilk + Em^fc Qjm 

< J2l^j,m^k llm + zZl^j Ilk + s j Hm^k Ijr 

Case (b), j = k: if j = k < p 

J2tyj,m^j llm +Sj E/^j Qlj + s j Hm^j Ijr, 

J2ljLj,m^j llm +Sj Ez^j Qlj + Em/j Qji 



Ez,m Hm 
+Sj ■ Sk ■ qjk 
+~Sk ■ Qjk 
+Sj ■ q jk 

Ez,m Qlm 
+Sj ■ Sj ■ qjj 

+ s i ' Qjj 



, J2l^j,m^j Qlm + Qlj 

Case (c): if j < p < k 



+ s j Em^j Qjm +Sj ■ qjj 



-I- [£] 

-i ■ Im J-l 

Liw J 



Ez^- «fc 

Em^fc 

. Ez,£j,m^fc 9Zm + Sj Em^fe 9j m 
Case (d): if k < p < j 

Qjk z 
J2l^j Qlk -- 
Ylrni^k Qjm z 
J2l^j,m^k Qlm + s k Ez^j Qlk z 



-- 

: 




-i • Im [4-1 . 

Liw J 



Proof. The sum ^"m=l QimW{F\ ' ) •( x )(^ 7 a '"^fc^') * s explicitely written in the different cases 
(a), (b), (c) and (d) as it was done in Lemma [4.61 Then the linear independence of the following 
families of functions is used to get the above systems for the gz m 's: 



THE MULTIPLE TUNNEL EFFECT ON A STAR-SHAPED NETWORK 15 

• Case (a): e Ax+Bx \ e Ax ~ Bx ' , e - Ax+Bx ' ', e - Ax ~ Bx> , 

• Case (b): cos(Ax) cos(Bx'), cos(Ax) sin(Bx'), sin(Ax) cos(Bx'), sm(Ax) sin(Bx'), 

• Cases (c) and (d): e Ax cos(Bx'), e Ax sm(Bx'), e~ Ax cos(Bx'), e~ Ax s'm(Bx'), 

where A and B are any fixed non- vanishing real numbers. 

On the other hand: if the qi m satisfy the system indicated in the lemma, on both sides of the 
equation of the lemma we have the same linear combination of the functions given above. Thus 
equality holds. □ 

Having the linear systems in the above lemma at hand, it remains to solve them. Indeed this 
is possible and we briefly indicate the necessary steps. 

Firstly, if A < ai, we have £j G iR for all j G {1, . . . , n}. Thus, for all j, k G {1, . . . , n} and all 

(x, x') G Nj x N k 



Im 



Now, let A be fixed in (a p , a p+ i), with p G {1, . . . , n}, remembering the convention a n+ \ = +oo. 
Due to the first equation of the corresponding system of cases (a), (c) and (d), the matrix 
Q '■= ( ( ?Zm)" m= i has to be of the form 



Qp 












with a p x p matrix Q p . 

Hence, the equations Yli^j Qik = and ^2 m ^ k Qjm = in case (a) are obviously fulfilled, as 
well as Ilk = in case (c) and Ylm^kQjm = in case (d). This means that only three 

equations remain from the cases (a), (c) and (d): 

£i, m ?im = Im[i] 
E m ^kHm + SkEQik = -Mm[i] 



Qjm 



Im [J-l , 



where in all the sums, I and m belong to {1, . . . ,p}. Now it is important to note that these 
three equations are already contained in the following system corresponding to the case (b), 
i.e. the only conditions for the g; m 's are the four following equations for a fixed (j, k) such that 
j < k < p or k < j < p: 



J2l^j, m ^k 9lm +$k Ilk +$j J2m^k 

J2l^j,m^k Qlm +Sk J2l^j Qlk + J2m^k Qjm 
+ Y^l^j Qlk 



- J2l^j,m^k Mm 

and, for j = k < p: 



~52l,m Qlm 

+Sj ■ Sfc • qjk 
+~Sk ■ Qjk 

+ S j Sm^fc Qjm + s j ' Qjk 



^2l,m Qlm 



J2l^j,m^j Qlm + s j J2l^j Qlj + s j Hm^j Qjm +Sj ■ Sj ■ qjj 
J2l^j,m^j Qlm +Sj Qlj + Ylm^j Qjm +Sj ■ qjj 

J2l^j,m^j Qlm + J2l^j Qlj + s j Ylm^j Qjm +$j ' Qjj 

where, once more, I and m belong to {1, . . . ,p} in all the sums. 



- Im [£] 



Ma 
- im ttf] 

— i ■ Im T J- 
i ■ Im [ J-l . 

iiw J 



(14) 



(15) 
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Now, we denote by Q jk := Qlk and Q' jk := Y. m ^k U™- Using the fact that Im [^] = 

Re , the above system (fT3|) is 

Yul^j^k Hm = Im [h] - Qjk - Q'jk - Ijk 

(sk-l)Qjk+ (sj ~ l )Q'jk+ {sjSk~-l)qjk = -2Re [i] 
(Sfc-l)Q ifc + (5fc-l)<&k = -ilm[i] -Re [i] = -i 

Since s,- — 1 = — the last three equations may be rewritten as 

(-icj£jw)Qjk+ {ickCkw)Q' jk + (-icj£jW + ick£kW -\w\ 2 )q jk = (2c^c k £j£ k )Re [±] 

Qjk+ Qjk = f^J 

Qjk + 9jk = j^p 

and the Gauss method yields that qj k must vanish for j ^ k. 

In the case j = k, we rewrite system f)15[) as above and the three equations to be solved turn 
out to be 



{-icj£jw)Qjj+ (icj£ j w)Q' jj + (-icjtjw + icjtjw - \w\ 2 ) qjj = (2c^) (Re [X] - ^- 
'•■• i.i ' '/./'./' ~~ ]™p 

The Gauss method once more gives the only possible solution qjj = for any j £ {1, . . . 

With this candidate for a solution at hand, it is pure calculation to show that it is indeed a 
solution. Thus we have shown the following result. 

Theorem 5.3. Let A be defined as in Section^ and the generalized eigenf unctions F^'° be given 
by Definition \3.4\ Take f £ H = YYj=i L 2 (Nj) , vanishing almost everywhere outside a compact 
set B C N, and let -co < a < b < +oo. Then for all x £ N 

(E(a,b)f)(x)= / fe V^(A)F A ~''(x) / F^(x')f(x')dx' dX, 



1=1 



where 



JO, if\<ai, 

llW '— \ c,(A)g,(A) 



— . if ai < A. 



Furthermore, for almost all A £ K the matrix <27 jm (A) — <^ m l( ajj+00 ) (A)q(A)6(A)/HA)| 2 * s the 
unique matrix satisfying (fT2]) . 



6. A DIRECT APPROACH 

We have seen in the preceding section that a matrix g(A) satisfying (|12p exists and is unique 
up to a null set. It is the aim of this section to deduce an alternative representation for q(X), 
involving only n x n-matrices and not an (3n 2 + 1) x n 2 system as above. Since this approach is 
essentially independent of the special setting, it should be more convenient for generalizations. 

In the following, let us consider the complex-valued generalized eigenfunctions F^' k as func- 
tions on N and not as elements of n?=i C°°(Nj). We introduce the notation 
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x 



Denoting by e k = (d~ik)i=i,...,n = (0, . . . , 0, 1, 0, . . . , 0) T the k-th unit vector in C™, we set di(X) := 
F\(Q) and for j = 2, . . . , n we fix Xj G Nj and set dj(X) := F\(xj). Due to the form of our 
generalized eigenfunctions, we then have 

n 

d 1 (X) = ^2e k = (l,...,l) T , 
k=i 

dj(X) = (3jej +otj^e k = (ay, . . . ,oy, f3j,atj, . . .,aj) T for j = 2,... ,n 

for suitable ay,/3j G C. 

Using these vectors, we now define 

n 

D(X) := J>(A)eJ. 

i=i 

Since ay ^ f3j, for every j G {2, . . . , n}, by construction, the matrix D is invertible for any choice 
of (x2, ■ ■ ■ , x n ) provided that ay ^ for all j G {2, . . . , n}. Indeed det D = \\^ =2 {Pj — ctj). 
Denoting by C(A) the diagonal matrix with c\\ = i and cyy = ictj for any j G {2, . . . , n}, we can 
formulate our theorem. 

Theorem 6.1. The matrix q(X) := {qi m W)l,m=i,...,n satisfying (fT2|) is given by 

1 n 

9(A) = (^A)" 1 ) 7 Im(— y ^^(Af^+i^CAjJ^W 1 ) 2, 



= ( J D(A)- 1 flm(-^yC(A) J D(A))p(A)" ) T 

for almost all X > a\. As previously, j has to be understood modulo n, that is to say, j + 1 = 1 
if j = n. 

Proof. By Lemma 15.11 the function q satisfies (|13p . Using the matrices and vectors introduced 
above, for fixed j G {1, . . . , n} this can be rewritten as 

Im(-^yF A (x) T ei+1 eJF A (x , )) = F x (x) T q(X) F\(x') 

for x G Nj and x' G N. 

Setting x = Xj and x' = x k we obtain by the definition of dj 

±-dj(X) T e j+1 e]d k (X)) = dj(X) T q(X)dkJX) 

for all j, k G {1, ... , n} and thus 

n . n ^ ran 



I,n U(A) 



k=ij=\ v 7 fc=lj=i 
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Using Ylk=i dk^k = we obtain 

Since D is invertible and the relation Y^ 1 j=i e jdj{X) T ej + ieJ = —iC{\) holds true, we now get 
the desired formula for q. □ 



7. A Plancherel-type Theorem 

In this section we prove a Plancherel type theorem for a Fourier type transformation V 
and its inverse V~ l associated with the generalized eigenfunctions (F7' k )k=i ... n introduced in 
Definition 13.41 Due to the fact that V and V~ x turn out to be not analogous, the surjectivity 
of V will not follow from its injectivity and needs a separate proof. 

Furthermore, we show that the fact that a function u belongs to the space D(A k ) can be 
formulated in terms of the decay rate of Vu. These results should be useful for the resolution of 
evolution problems involving the spatial operator A as well as for the analysis of the properties 
of a solution. For this part, we follow Section 4 of [2]. 

Definition 7.1. Let qk be the coefficients defined in Theorem 15.31 Considering for every k = 
l,...,n the weighted L 2 space L 2 ((ak, +oo), qk), we set L 2 := YYk=i L 2 {{ak, +oo), qk). The 
corresponding scalar product is denoted by 

n „ 

(F, G) q -=J2 qk(X)F k (X)GkJ\) dX 

k= \ J (Ofc:+0°) 

1/2 

and the norm by |-F| 9 := (F,F) q . 

Note that L 2 , is a Hilbert space, since it is the product of the Hilbert spaces L 2 ((ak, +oo),qk). 

Now, we define the transformation V, together with its inverse Z, which is later on shown to 
diagonalize A. 

Definition 7.2. i) For all / £ L l (N, C) and k = 1, . . . ,n the function (Vf)k : [o*., +oo) -> 
C is defined by 

(Vf)k(X) := I f(x)(F-> k )(x) dx. 
Jn 

ii) We choose a cut-off function % G C°°(M) such that x = on (— cxD,a n + 1) and % = 1 
on (a n + 2, +oo). Then we set X to be the space of all (G?i, . . . , G n ) G L 2 , such that 
Gk G C°°([afc, +oo), C) and xGk can be extended by zero to an element of <S(R) for all 
k G {l,...,n}. 

Now, for every (G\, G n ) G X we define Z(G?i, . . . , G n ) : iV C by 
Z(G 1: G n )(x) :=J2f q k (X)Gk(X)(F-' k )(x) dX. 

k=1 J (a fc ,+oo) 

Lemma 7.3. Let f G L 2 (N) and G = {G x , . . . , G n ) G X. T/ien G G L 2 , Z(G) G L 2 (iV) = H, 
Vf = ((Vf) 1} . . . , (Vf) n ) G L\, and 

(G,Vf) q = (Z(G),f) H . 
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Proof. We have 

n 

L 2 ((afc, +00), gj.) 



fc=l 

n 



E 



(a fe ,+oo) 



q k (X)G k (X)(Vf) k (X) dX 



V / q k (X)G k (X)([ f(x)(F-> k )(x) dx) dX 

^J(a fc ,+oo) VAT ' 

n „ 

Y, / q k (X)G k (X){F-> k )(x)dx)W)dx = (Z(G)J) H . 



k =i 



(a k ,+co) 



(16) 



It only remains to make sure that the assumptions for Fubini's Theorem are satisfied to justify 
(fT6j) , In fact, it is sufficient to estimate (X,x) \-t q k (X)(F^ ,k )(x) on (a k ,a k+ i) x N for a fixed 
fc € {1, . . . , n}, since / is compactly supported and G k is rapidly decreasing. 

In order to do so, recall that c k ^ k (X) = \J c k (X — a k ) belongs to 1R if and only if A > a k and 
to iM. otherwise and that for any A > a%, we have |u>(A)| 2 > XT=i c i\^ ~ a i\ due t° Lemma [3771 
Thus, putting in the expression for s k (X), cf. Definition 13.41 we get for x G N k 



\q k (X)(F^ k )(x)\ 



c k £ k {X) 



cos 



< 



HA)| 2 

\c k £ k (X)\ 



^k(X)x) - is k (X)sm(£ k (X)x)) 



< 



\ckCk(X)\ 



Er=i c /i A - oil 



Ya=iCi\X- a t \ |c fc £ fc (A)| 
Furthermore, for x G Nj, j > k, we have 



E?=iciMX)\ < E?=iV^iV\^\ 



H?=iCi\X-ai\ 



Cl + |s fc (A)|) 



(17) 



\q k (X)(F-' k )(x)\ 



Finally, for x G Nj, j < k, 



\w(X)\< 



< 



exp 



I dj — X 



\q k (X)(F-> k )(x)\ 



\w(X)\ 2 



< 



Ck\/\X - a k \ 



(18) 



(19) 



For all three cases the bound is an L\ function of A, which is enough for the application of 
Fubini's Theorem, due to the properties of / and G. □ 

Lemma 7.4. Consider f G flLi C ™( N k)- Then {{V f\, {V f) n ) G X. 

Proof. As in Section 4 of [2], (Vf) k is a linear combination of Fourier and Laplace transforms 
of functions in C^(Nj), j G {1, . . . , n}. □ 

The next step is to show that V is an isometry and Z is its right inverse. The proof of the 
following theorem is analogous to that of [2j Theorem 4.9]. 



Theorem 7.5. Endow Ylk =1 C™(N k ) with the norm of H = \\l =l L 2 {N k ). Then 

i) y : 11^=1^0° C^fc) ~~ ^ L 2 , is isometric and can therefore be extended to an isometry 
V:H^L 2 . 

In particular, for all f G H we have \Vf\g = (f, 

ii) Z = V- 1 on V(H). 

iii) Z can be extended to a continuous operator on L 2 
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In the sequel we shall again write V for V for simplicity. Note that thanks to the density of 
X in Lq, from Lemma 17.31 we get the Plancherel type formula 

(G, Vf) q = (Z(G), f) H for all / G L 2 (N) and G G h\. (20) 

Set 

Y := {(Gi, . . . , G n ) G L 2 q : supp(Gj) compact and supp(Gj) n {ax, . . . , a n } = 0, j = 1, . . . , n}. 

Then Y is a dense subspace of and for our subsequent considerations it is important to know 
that also for all G G Y the transformation Z is given by the formula from Definition 17.21 ii) . 
This is the aim of the following lemma. 

Lemma 7.6. Let G G F. Then 

Z{G)(x) = W q k (\)G k (X){F-> k )(x) dX. 

k=1 J(a k ,+co) 

Proof. Let / G Uk=i C Z°{ N k)- Then, using ([20]), we get 

(Z(G),f) H = (G, Vf) q = J2f G k (X) f F-' k (x)f(x) dx q k (X) dX. 



In order to apply Fubini's Theorem, we estimate for every k = 1, . . . , n 

G k (X)F-> k (x)f(x)q k (X)\dx=\G k (X)q k (X)\ [ \F-' k (x)\\f(x)\ dx 
N JN 

< \G k (X)\ sup \qk(X)F-> k (x)\-\\f\\ LHN) . 

xGsupp(/) 

Investing the estimates from (I17p - (|19p . we can estimate this further by 

< c J: — \G k (X)\ -\\f\\ LHN) , 

Vl A - a k\ 

and since G G Y this function is in L l ((a k , +oo)), so we may apply the Fubini Theorem and 
obtain 

(Z(G),f) H = f J2 f q k (X)G k (X)F-' k (x) dX f{x) dx 

JN k=1 J[a k ,+oo) 



ex 



k= l J(a k ,+co) 



q k (X)G k (X)F-' k dXJ 



H 



This implies the assertion. □ 

We now want to show that V diagonalizes the operator, i.e. V(A J u)(X) = X J V(u)(X) for all 
u G D(A J ). In order to formulate this precisely for a measurable function ip : R — > R we denote 
the corresponding multiplication operator on l? q by M^, i.e. for a function F G l? q we have 

{M^F) k {X) = tp(X)F k (X), AG[a fe ,+oo), k G {1, . . . , n}. 

Then we have the following lemma, whose proof is again analogous to that of [21 Lemma 4.12]. 

Lemma 7.7. i) Let j G N and pj : R — > R be defined by Pj(x) = x 3 . Then for any 
f G D(A>) we have 

V(Aif) = M Pj Vf. 
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ii) Let ip : [ai,+oo) — > R be a bounded, measurable function. Then ip(A) defined by the 
spectral Theorem is a bounded operator on H and for all f G H we have 

V^(A)f) = M^Vf). 

In order to show injectivity of Z, and thus the surjectivity of V, we will need the following 
property of the kernel. 

Lemma 7.8. Let G G ker(Z) and let ip : [ai, +00) be measurable and bounded. Then M^G G 
ker(Z). 

Proof. Due to the boundedness of t/j, it holds D(ip(A)) = H and M^G G L? q for G G L? q . There- 
fore, if / G L 2 (N) with compact support and G G ker(Z), we have using (|20|) and Lemma 17.71 

(Z(M^G)J) H = (M^G,Vf) L i = {G,Mi,Vf)q = (G, Vip(A)f) L 2 
= {Z{G)MA)f) H = (0, 4>(A)f) H = . 

□ 

This lemma allows to reduce the proof of the injectivity for Z to the frequency bands (a^, a^+i) 
for k = 1, . . . ,n (set again a^+i = +00). There the injectivity of the Fourier sine transform will 
be used. 

Theorem 7.9. ker(Z) = {0} and thus V is a spectral representation of H relative to A. 

Proof. Let G G L? with Z{G) = be given. We fix j G {1,... , n} and consider a compact 
interval I C (a,-, Oj+i). Since ljG lies in Y, invoking Lemmas 17.81 and 17.61 we obtain for almost 
all x G N 

n „ j ■ 

= £(1/G)(x) = / % (A)l 7 (A)G fc (A)F-' fc (x) dA = £ / % (A)G fc (A)F"' fc (x) dA. (21) 

fc=1 J(a fc ,+oo) fc=l^ 

In this equation, we now let x to zero. This is possible, since the hypotheses of Lebesgue's 
Theorem are satisfied thanks to the presence of 1/. This implies thanks to F7 ,k (0) = 1 



°= / y>fc(A)G fc (A) dA. 



Since this is true for every compact interval / C (a 3 -,aj+i) and since g^G^ G L\ QC ((aj, Oj+i)) for 
all A; = 1, ... , j, we can argue as in the proof of Lemma 15. II to obtain 

3 

%(A)G fc (A) = for a.a. A G (aj,a j+1 ). (22) 

k=l 

In a second step, we now consider x G N% for some £ < j in (|2ip . This leads to 

= ^^g fc (A)G fc (A)exp(-^(A)x) dA + Jto{X)G t (\) [cos(&(A)z) - i«j sin(&(A)x)] dA 



= 12 / %(A)G fc (A)exp(-i&(A)x) dA - i / q e (X)G e (X){s e (X) - l) sin(&(A)x) dA. 

Thanks to ()22[) . the sum in the last expression vanishes, so we are left with 

= / 5<(A)G*(A)(^(A) - 1) sin(&(A)x) dA. 
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In this case &(A) = y(A — ae)/ce is real, so we may substitute p = &(A) and obtain 
= 2q / <^(a,£ + c e p 2 )Gt(ai + cip 2 ) [se(a,£ + Qp 2 ) - l] sin(px)p dp 

= 2q / (p)^^ 1 ^))^^" 1 ^)) kfe _1 0)) " l]psin(px) dp. 

J(0,+oo) 

This last integral is now the Fourier sine transform of l^m (qe C^X^ ° ^i" 1 ) [( s ^ ° — 
so by the injectivity of this transformation we find 

2Q^(^ 1 (P))G'K^ 1 (P))[(S£(^ 1 (P)) - l]p = for a.a. p G 

Dividing by p and resubstituting this implies 

2c*<fc(A)G<(A)|>*(A) - 1]&(A) = for a.a. A 6 J. 

Observing, that and — 1 all may only vanish in the points ai, . . . , a n , this finally implies 
G^(A) = for all £ G {1, . . . , j} and almost all A G (a-,-, Oj+i). Since j G {1, . . . , n} was chosen 
arbitrarily, the assertion follows. □ 

Note that these considerations also provide a proof for the special case n = 2, whose proof in 
[3 Theorem 4.17] is not entirely correct. 

Finally, we get the following characterization of the spaces D{A ) ). 

Theorem 7.10. For j G N the following statements are equivalent: 

i) u G D(AP), 

ii) A i — y \i(Vu){\) G L\, 

iii) A h-> A J (yw)fc(A) G L 2 ((afc, +oo), / or all k = 1, . . . ,n. 

Proof. Due to Theorem 17.51 it holds (A>u, A^u)h = IV^-A^u)^. Now Lemma 17.71 implies 

□ 

The above results provide explicit solution formulae for evolution equations involving the 
operator A. For example for the wave equation ii + Au = with u(0) = uq and u(0) = 
a formal solution is given by u(t) = Z cos(\^Xt)VuQ. Our aim in the future will be to study 
properties of these solutions, as indicated in the introduction. 



8. Appendix 

The book of J. Weidmann |27] describes a general approach to the spectral theory of systems of 
Sturm-Liouville equations, which is in principle applicable to our setting. In this point of view, 
our problem is seen as a system of n equations on (0, +oo) coupled by boundary conditions in 
0. Thus the kernel of the resolvent is a matrix-valued function /C(-, -, A) : (0, +oo) x (0, +oo) — > 
C nxn . The relation to Definition 13.61 is given by 

ICij(x,x',X) = K(x,x',X) where (x,x') G 2Vj x Nj = (0, +oo) x (0,+oo). 
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The fundamental hypothesis of Weidmann is that the kernel of the resolvent can be written in 
the following form: 



fC(x, x , X) 



p 

Yfq=l (S2 aic > Wl ^ Xj A )) W <?( X '> A ) T ' fOT X> - X ' 
1=1 



=: m?(x,\) 

P q (23) 

Eg=i (X)A<2 wi(x,X)) w q (x',X) T , for x' > x, 
1=1 



= : mq(x,X) 

where ai q ,(3i q G C and the w q : [0, oo) x C — > C n are such that {w q (-,X) : q = l,...,p} 
is a fundamental system of ~ker(Af — XI), the space of generalized eigenfunctions of A. Here 
Af : D(Af) — > C (M) is the formal operator, in our case Af = A = (— • d% + ak)k=i,...,n but 
D(Aj) = YYk=i C 2 (-^V"fc)) i.e. the operator A without transmission conditions nor integrability 
conditions at oo. Clearly in our case dim(ker(Aj — XI)) = 2n and thus p = 2n. 

In contrast to the typical applications treated in [27] as for instance the Dirac system, we 
consider in this paper a transmission problem, i.e. intuitively the components of the w q are 
functions on different domains (while mathematically all are equivalent to (0,+oo)): 
the branches of a star. For all such applications it is highly important to use only generalized 
eigenfunctions satisfying the transmission conditions (To) and (T\), for example Theorem 17.101 
would be impossible otherwise. 

Supposing the ansatz (|23p of Weidmann, this is not possible, what we shall show now. 

In Definition 13. 61 we have given an explicit expression for the (unique) kernel K of the resolvent 
of A, using only elements of ker(Ar — XI), where At : D(At) — > H satisfies At = A = {—Ck-d^.+ 
Ofc)fc=l,...,n and D{A T ) = Uk=i C 2 (N k )n{{u k )l =1 satisfies (T ), (Ti)}. Note that ker(^ T - XI) is 
the n-dimensional space of generalized eigenfunctions of A satisfying the transmission conditions 
(To) and (Ti), but without integrability conditions at infinity. 

Suppose that we have a representation of K as given in (j23|) satisfying 

w q (;X) £ker(A T -XI), q = l,...,p, X G p(A) 

and thus we can take p = n. Let us fix x G iVi and A G p(A). Then the m q (x, A) =: m q G C are 
constants and the expression 

p 

g{x',X) :=^2m^w q (x',X) 

9=1 

defines a function 

g(-,X) : {x' G [0,oo) : x' > x} -> C n . 

Clearly, the functions 

[w q (;X)]j :Nj^C 

can be uniquely extended to entire functions in x' because they are linear combinations of 
e ± l 0( A ) x . Thus the [g(-, X)]j are entire. Our hypothesis w q (-,X) G kei(AT — XI), q = 1,... ,p, 
implies thus 

g(-,X) eker (A T -XI), q=l,...,p. (24) 
But comparing (|23p with Definition 13.61 yields 

g{x', A) = F^ 2 (x) ■ F}'\x') for x' G JV 2) . . . , 
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and 

g{x', A) = F^' 1 (x) ■ F±' 2 (x') for x' G N x , 

which makes sense after analytic continuation. Using the assumption that g satisfies (To), we 
get from these two equalities, putting x' = £ N2 into the first and x' = G N% into the second 

F±> 2 (x) = F^(x)F^(0) = 5 (0, A) = F^(x)F^ 2 (0) = F^(x). 

But inspecting the definitions of the generalized eigenfunctions, bearing in mind that x E N\ 
was arbitrary, this implies 

cos(£i(A)x) ± zsin(£i(A)x) = cos(£i(A)x) ± isi(X) sin(^i(A)x) 

and thus si(A) = 1. This finally implies Y^j=i c j£jW = 0> which is impossible for all real A < a%, 
since then £j(A) is purely imaginary. Furthermore, Y^=i c j£j * s analytic, so this equality can, if 
ever, only be fulfilled on a discrete set of A G C. 
We have thus proven: 

Theorem 8.1. Representation (|23p of the kernel JC of the resolvent of A is not possible, using 
exclusively generalized eigenfunctions 

w q (-, A) G ker(A T -XI), q = l,...,p. 

The reason for this is, shortly speaking, the rigidity of (f23j) caused by the use of the same 
linear combinations of generalized eigenfunctions above the diagonals of all Nj x iV& (idem 
below). This is not compatible with the fact that the kernel is non-smooth only on the main 
diagonals, i.e. the diagonals of Nj x Nj, cf. Definition 13.61 and Figure [TJ 

This means that the approach of J. Weidmann in [27], when applied to problems on the star- 
shaped domain, does not sufficiently take into account its non-manifold character. The resulting 
expansion formulae would use generalized eigenfunctions which are in a sense incompatible with 
the geometry of the domain. This would have undesirable consequences: for example, the 
important feature of Theorem IT. 101 that the belonging of u to D(A J ) is expressed by the decay 
of the components of Vu would be impossible, due to artificial singularities in the expansion 
formula. 
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